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Abstract. Day and night airborne thermal infrared image data at 5 m spatial 
resolution acquired with the 15-channel (0-45 /im- 12*2 /im) Advanced Thermal 
and Land Applications Sensor (ATLAS) over Alabama, Huntsville on 
7 September, 1994 were used to study changes in the thermal signatures of urban 
land cover types between day and night. Thermal channel number 13 
( 9-60 /mn-1 0*2 /an) data with the best noise-equivalent temperature change 
(NEAT) of 0-25 C after atmospheric corrections and temperature calibration 
were selected for use in this analysis. This research also examined the relation 
between land cover irradiance and vegetation amount, using the Normalized 
Difference Vegetation Index (NDVI), obtained by ratioing the difference and the 
sum of the red (channel number 3: 0-60-0-63 /mi) and reflected infrared (channel 
number 6: 0-76-0-90 /im) ATLAS data. Based on the mean radiance values, stand- 
ard deviations, and NDVI extracted from 351 pairs of polygons of day and night 
channel number 13 images for the city of Huntsville, a spatial model of warming 
and cooling characteristics of commercial, residential, agricultural, vegetation, 
and water features was developed using a GIS approach. There is a strong 
negative correlation between NDVI and irradiance of residential, agricultural, 
and vacant/transitional land cover types, indicating that the irradiance of a land 
cover type is greatly influenced by the amount of vegetation present. The predom- 
inance of forests, agricultural, and residential uses associated with varying degrees 
of tree cover showed great contrasts with commercial and services land cover 
types in the centre of the city, and favours the development of urban heat islands. 
The high-resolution thermal infrared images match the complexity of the urban 
environment, and are capable of characterizing accurately the urban land cover 
types for the spatial modeling of the urban heat island effect using a GIS approach. 


1. Introduction 

The main objective of this research is to apply airborne high-resolution thermal 
infrared imagery to study the phenomenon of the urban heat island, using Huntsville, 
Alabama, a medium-sized American city, as the study area. The occurrence of urban 
heat islands represents human-induced urban/rural contrast, which is caused by 
deforestation and the replacement of the land surface by non-evaporating and non- 
porous materials such as asphalt and concrete. The result is reduced evapotranspir- 
ation and more rapid runoff of rain water (Carlson 1986, Kim 1992). The urban 
landscape forms a canopy acting as a transitional zone between the atmosphere and 
the land surface. The composition and structure of this canopy have a significant 
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impact on the thermal behaviour of the urban environment (Goward 1981 ). In other 
words, apart from the properties of urban construction materials, the layout of 
buildings and pavements within an urban landscape can also affect the appearance 
of the urban heat island (Oke 1987). 

Research on the trends of surface temperatures at rapidly growing urban sites in 
the U.S.A. during the last 30 to 50 years suggests that significant urban heat island 
effects have caused the temperatures at these sites to rise by 1 to 2°C (Cayan and 
Douglas 1984, Karl et al. 1988). Urban heat islands have caused changes in urban 
precipitation and temperature that are at least similar to, if not greater than, those 
predicted to develop over the next 100 years by global change models (Changnon 
1992). 

Satellite remote sensing, particularly NOAA AVHRR data, has been used in the 
study of urban heat islands (e.g., Lee 1993, Matson et al. 1978, Brest 1987, Gallo 
et al. 1993, Roth et al 1989). Because of the low spatial resolution (M km at nadir) 
of the AVHRR data, these studies can only examine and map the phenomenon at 
the macro-level. As Roth et al. (1989) observed, the heterogeneity of the urban 
environment rendered the definition of this surface very difficult, and the smoothing 
of signals within a pixel introduced bias in the radiant temperatures extracted. 
Recently, Nichol (1994) has successfully utilized Landsat Thematic Mapper (TM) 
thermal infrared data (with a spatial resolution of 120m) to derive surface temperature 
data for some housing estates in Singapore. But the spatial resolution of the TM 
thermal infrared data is still inadequate to capture all the complex temperature 
changes of the urban environment. 

This research utilized 5 m thermal infrared data acquired specially from an aircraft 
to characterize more accurately the thermal responses of different land cover types 
in the urban environment as input to urban heat island detection. 

2. Data acquisition and processing 

The thermal infrared data were acquired by the Lockheed Engineering and 
Science Company as Mission M424 over the greater Huntsville, Alabama area, on 
7, September 1994, a clear day with less than 5 per cent of cloud cover in the sky, 
using the Advanced Thermal and Land Applications Sensor (ATLAS) sensor system 
on board a NASA Stennis LearJet. The ATLAS is a 15-channel imaging system 
which incorporates the bandwidths of the Landsat Thematic Mapper with additional 
bands in the middle reflective infrared and thermal infrared range (table 1). All 
15-channel data were acquired on the same date at 10 m and 5 m resolutions. To 
optimize the detection of warming and cooling of urban land surfaces, the image 
data were planned to be acquired first around solar noon (12:00-13:00) and then 
repeated 2-3 h after sunset. In all, starting from 11:00 Central Time, nine daytime 
flight lines were flown from an altitude of 5000m to image at 10m resolution, and 
six daytime flight lines were flown from an altitude of 2500 m to image at 5 m 
resolution. All these flight lines were repeated at night beginning at 20:30. Central 
Time. Colour infrared aerial (CIR) photography was also acquired simultaneously 
using a Zeiss aerial camera (with 152 mm focal length and 23 cm x 23 cm photograph 
format) during the daytime flights. Ground data collection, including the establish- 
ment of GPS ground control points, measurement of surface temperatures at selected 
sites of different cover types throughout the city, and the concurrent launches of 
radiosonde were also carried out at the time of overflights. In addition, meteorological 
station data were obtained from the Army at Redstone Arsenal and at Marshall 
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Band width NER NEAT MTF @ 


Channel 

limits (/fm) 

mW/cm 2 (/im) 

°C 

2 mrad 

Cooling 

1 

0-45-0*52 

<0-008 

N/A 

0*5 

Ambient 

2 

0*52-0-60 

<0-004 

N/A 

0-5 

Ambient 

3 

0-60 0*63 

<0*006 

N/A 

0-5 

Ambient 

4 

0-63-0-69 

<0*004 

N/A 

0-5 

Ambient 

5 

0-69 0*76 

<0-004 

N/A 

0-5 

Ambient 

6 

0-76-0*90 

<0-005 

N/A 

0*5 

Ambient 

7 

1*55-1*75 

<0-05 

N/A 

0-5 

77 K 

8 

2-08-2*35 

<0-05 

N/A 

0-5 

77 K 

9 

3*35-4-20 

N/A 

<0-3 

0-5 

77 K 

10 

8*20-8-60 

N/A 

<0*2 

0-5 

77 K 

11 

8*60-9-00 

N/A 

<0*2 

0*5 

77 K 

12 

900-9-40 

N/A 

<0-2 

0-5 

77 K 

13 

9*60-10*2 

N/A 

<02 

0-5 

77 K 

14 

10*2-11-2 

N/A 

<0*2 

0-5 

77 K 

15 

11*2-12*2 

N/A 

<0*3 

0-5 

77 K 


Space Flight Center during overflights. The Portable Ground Atmospheric 
Measurement System (PGAMS) was also used to collect data for atmospheric 
correction. Altogether over 10 Gb of digital image data were obtained by this flight 
mission. Despite some glitches in channels 7, 8, and 9 (the mid-infrared channels) 
and some smears in channel 3, the overall quality of these digital data, particularly 
the thermal infrared data of channels 10 to 15, was judged to be very good. The 
CIR photography is excellent in quality, and, because of its superior spatial reso- 
lution, provides ‘ground truth’ information against which the digital image data can 
be checked. 

For this particular research, the focus has been on the six thermal infrared bands 
(channels 10 to 15, in the spectral range from 8-20 /rm to 12-2 fim; table 1). These 
data have recorded long-wave radiation emitted by both the natural and artificial 
surfaces in the city of Fluntsville. The original image data recorded the radiation as 
digital number (DN) in 8-bit format with integer values ranging from 0 to 255. 
Before these image data can be used, they have first to be corrected for the attenuation 
effect of the atmosphere, namely, transmittance and path radiance, and second, to 
be calibrated to produce accurate temperature measurements. Atmospheric correc- 
tion was achieved by applying the MODTRAN program developed by the United 
States Air Force Geophysics Laboratory, which estimates atmospheric transmittance 
and radiance for a given atmospheric path at moderate spectral resolution over the 
operational wavelength region from 0-25 to 28-5 /tm calculated in lcm -1 frequency 
interval (Berk et al. 1989). The input parameters to the program for the atmospheric 
correction are the radiosonde data of atmospheric profiles during the ATLAS over- 
flights. Temperature calibration of the ATLAS sensor for each channel was achieved 
by using onboard low and high temperature blackbodies which were referenced at 
the beginning and end of each scan line. With a knowledge of the emittance value 
for the blackbodies, Planck’s equation was applied to calculate the ATLAS reference 
blackbody radiance for each channel of the sensor system. By combining the outputs 
of the MODTRAN atmospheric correction program with the high and low black- 
bodies temperature calibration of the ATLAS sensor system, the Earth Resources 
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Laboratory Applications Software (ELAS) module TRADE (TIMS < CAMS, 
ATLAS > Radiant Energy) computed the system transfer equation (STE), within- 
channel normalization, atmospheric correction, and the conversion of the 8-bit DN 
of each pixel of the image data into 32-bit atmospheric radiance in units ofWcm 2 
sr 1 /an 1 (Graham et al. 1986, Luvall et al. 1990, Anderson 1992). 

In this paper, to provide a better appreciation of the thermal characteristics of 
the land cover type, radiance is converted into irradiance, or radiant flux density 
incident on a surface in Wm“ 2 (Monteith 1973). The use of energy values avoids 
the need for emissivity correction for land cover in the case where temperature values 
are used. 

This research will make use of the high-resolution 5 m ATLAS data only because 
the complexity of human activities in the urban environment requires the highest 
possible spatial resolution for an accurate characterization of the urban land cover’s 
thermal responses. In the next phase of the research, the 10 m resolution data will 
also be employed for comparative study and to determine the optimum spatial 
resolution. In view of the need to produce a false colour composite image and to 
compute the Normalized Difference Vegetation Index (NDVI), four channels of data: 
channel 2 ( 0*52-0-60 jum: visible green), channel 3 (0-60-0*63 ^m: visible red), chan- 
nel 6 (0-76-0-93 /m: reflected infrared), and channel 13 (9-60-10-2 /un: thermal 
infrared) were extracted from the six flight lines acquired in the daytime, and one 
thermal channel (13) from the corresponding six flight lines acquired at night. 

Channel 13 (9-60-10-2 /im) is the thermal band of choice mainly because after 
the atmospheric correction and temperature calibration, channel 13 data exhibit the 
best noise-equivalent temperature change (NEzfT) of 0*25 °C (i.e., the temperature 
change across the target that would produce a signal-to-noise ratio of unity in the 
detector output). According to Wien’s Law, the surface temperatures recordable 
within the waveband limits of channel 13 are from 284 to 301 K. In other words, 
the waveband limits of channel 13 display the maximum energy per unit wavelength 
of terrestrial radiation. Image data from the other three channels: channel 2 
(0-52-0-60 jum), the green band, channel 3 (0-60-0-63 //m) the red band, and channel 6 
(0-76-0-90 /^m), the reflected infrared band, when displayed through the blue, green, 
and red guns of the computer monitor respectively as overlays, will produce a false 
colour composite image (simulating the colour infrared) which emphasizes vegetation 
vigour. When channel 2 data are displayed in green, channel 6 data are displayed 
in red, and channel 13, the thermal band, data are displayed in blue, the false colour 
composite image produced will also pinpoint hot and cool objects distinctly, in 
addition to indicating the vigor of the vegetation (figure 1). As for the night-time 
imagery, images can be recorded clearly only in the six thermal bands. A night-time 
image of channel 13 is shown in figure 2. 

3. The study area 

The city of Huntsville, located in north Alabama, is not a typical American city. 
Two dominating land uses within the city are the Redstone Arsenal and the Marshall 
Space Flight Center of NASA, both related to the U.S. Government. It is also the 
site for the campus of the University of Alabama in Huntsville (UAH). Services and 
industries are attracted here to serve them. Many of these are high-tech in nature, 
such as Intergraph, and their locations in the city are quite distinct. The city has 
developed spacious sites to locate these offices and industries in the form of a 
’research park’, developed over agricultural land. These offices are often intermingled 
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Figure 1. A false colour composite of the Madison Square Mall area of Huntsville, Alabama, 
formed by displaying channel 2 (0-52-0*60 /<m) image in green, channel 6 (0-76-0*90 //m) 
image in red, and channel 13 (9*60— 10*2 /an) image in blue. Colour changes from white, 
yellow, pink, orange, red, light blue, dark blue, and black, indicating changes of 
irradiance from very high to very low. 


with croplands of soybeans, cotton, and corn. The city does not have a large central 
business district (CBD). However, the street pattern in the CBD is more dense and 
distinct. The city is well served by highways and roads, the most important of which 
are Interstate 565, Highway 72, and Memorial Parkway (Highway 231/431). 
Topographically, the city is also quite unique. Located at an average elevation of 
236 m, the city’s eastern extent is limited by a north-south running barrier of 
mountains: Huntsville Mountain, Green Mountain, and Wallace Mountain. Its 
northern extent is marked by Smither Mountain and Wade Mountain. Southeast of 
the city is the Ward Mountain. Its southern extent is marked by the Tennessee River 
and its associated lakes and tributaries, where wetlands occur in abundance. It is 
interesting to know how such a combination of landforms and human activities has 
affected the development of the urban heat island. 

4. Analysis 

In order to make the best use of the high-resolution thermal infrared data, the 
characteristics of the thermal signatures of each land cover type in the city of 
Huntsville were studied. With the aid of the ELAS image processing software running 
on a Silicon Graphics (SG) workstation, the daytime and night-time images of the 
city was displayed simultaneously. A key of land cover types with codes was 
developed (table 2). From the displayed images, representative polygons of each land 
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Figure 2. Night-time channel number 13 (9*60-10-2 /^m) image of the Madison Square Mall 
area of Huntsville. Black indicates low irradiance and white indicates high irradiance. 

cover type were manually delineated from both the daytime and night-time images, 
in the same manner as training sets are delineated in a supervised image classification. 
The mean irradiance values and standard deviations of these polygons were extracted 
from the images. A total of 351 pairs of polygons were extracted. 

Because the urban heat island is related to the reduction of evapotranspiration 
from surface vegetation cover, it is useful to know the relationship between surface 
vegetation cover and water availability. Recent research has shown that NDVI is a 
good indicator of surface radiant temperature (Nemani and Running 1989, Gillies 
and Carlson 1995). NDVI for the city of Huntsville is therefore computed from the 
daytime image data using the following formula: 

NDVI = (ch6— ch3)/(ch6 + ch3) ( 1 ) 

where ch3 is the red band (0-60-063 fim) and ch6 is the reflected infrared band 
(0-76-0-90 fim) of the image data. The NDVIs were computed for each of the land 
cover polygons extracted from the daytime images. 

5. Results 

5.1. Day and night contrasts in irradiance 

The average irradiance values by land cover types are summarized in table 3. It 
is clear that during the daytime, commercial land cover types exhibit the highest 
irradiance (above 70 Wm“ 2 ), followed by services (65-70Wm -2 ), transportation 
(57-74Wm~ 2 ), and industrial uses (63-69 Wm’ 2 ). The standard deviations (SDs) 
of the daytime irradiance values are high, indicating the heterogeneous nature of 
these land cover types. By contrast, the lowest daytime irradiance is found in water 
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1. Residential (RE) 

a single family homes with >50% and of tree canopy 
b single family homes with between 50% to 10% tree canopy 
c single family homes with <10% tree canopy 
d multi-unit dwellings (apartments/condominiums) (generally 
without much vegetation canopy) 

2. Commercial (CO) 

a retailing (large) 
b retailing (small) 

3. Services (SE) 

a offices (high-rises): high density 
b offices (low buildings): low density 
c schools and universities 

d other institutions (government or private), e.g. NASA 

4. Industrial (MN) 

a light 
b heavy 
c warehouse 

5. Agriculture (AG) 

a cropland with row crops growing 
b cropland harvested/fallowed; not currently used 
c pasture 

6. Transportation (TR) 

a interstate highways 
b state highways 
c bridges 

d air fields (civilian and military) 
e parking lots 

7. Vegetation (VE) 

a coniferous forests 
b deciduous forests 
c mixed coniferous and deciduous 
d grass 

e wetland vegetation (visible) 

8. Water bodies (MA) 

a rivers and streams 
b ponds, lakes, and reservoirs 
c wetland (no vegetation visible) 

9. Recreational/open space (RC) 

a golf courses 
b parks 

c baseball fields/stadiums 

10. Vacant or transitional (VA) 


bodies (54-59 Wm~ 2 ) followed by vegetation (55-62 Wm 2 ) and agriculture 
(55-62 Wm -2 ). Residential uses occupy an intermediate position in irradiance 
(62-69 Wm^ 2 ) because they are composed of buildings and tree cover in varying 
degrees of mixture. 

At night, commercial, services, industrial, and transportation land cover types 
cool down rapidly to very similar irradiance values (45-47 Wm 2 ) so that they are 
not too much higher than those for vegetation and agriculture (table 3). Water is 
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Table 3. Average irradiance by land cover types for day and night. 


Land cover type 
(see key in table 2) 

Day 

Wm’ 2 

Day:SD 
Wm 2 

Night 

Wm‘ 2 

Night:SD 

Wm’ 2 

DifT. between day & 
night 
Wm’ 2 

1. Residential (RE) 

65-44 

8-21 

45*53 

1*97 

19*92 

REa 

61-85 

7*88 

45-13 

1-48 

16-72 

REb 

64-50 

8-78 

45-94 

1-92 

18-56 

REc 

65*98 

7*55 

44-55 

2-18 

21-42 

REd 

68-94 

8-33 

46-01 

2-30 

22-92 

2. Commercial (CO) 

70-90 

7*99 

46-26 

4-31 

24*64 

COa 

71-48 

7-56 

46-51 

4-33 

24-97 

COb 

70*37 

8*38 

46-03 

4*29 

24-33 

3. Services (SE) 

68-35 

6-99 

45*95 

3-40 

22-40 

SEa 

65-32 

8 13 

45-41 

1-97 

1991 

SEb 

69-27 

7-85 

44-92 

4-05 

2436 

SEc 

70-31 

7-85 

46-02 

3*30 

24-29 

SEd 

66*79 

6-24 

45*89 

3-43 

20-90 

4. Industrial (MN ) 

66-09 

6-55 

44-85 

3*28 

21-24 

MNa 

68*90 

7 10 

45-04 

3*57 

23-85 

MNb 

63*32 

4-56 

47-51 

1-63 

15-81 

MNc 

63-65 

6-90 

43 06 

3-80 

20-60 

5. Agriculture (AG) 

58-26 

I 09 

43-20 

0-51 

15-06 

AGa 

55*20 

1 10 

42-33 

0 51 

12-86 

AGb 

61 77 

1 00 

43-99 

041 

17-77 

AGc 

57-93 

1 14 

43-22 

056 

14-70 

6. Transport (TR) 

67-73 

4-53 

47*45 

1-69 

20*27 

TRa 

67*92 

7*52 

45*56 

3 19 

22-37 

TRb 

69*30 

3 96 

48-08 

1-38 

21-22 

TRc 

57-16 

2-75 

46-48 

204 

10-68 

TRd 

64-55 

5*23 

44-90 

2-43 

19-66 

TRe 

73-69 

4-41 

49-33 

1 05 

24-37 

7. Vegetation ( VE) 

56-90 

0-88 

45-56 

0-57 

10*34 

VEa 

54-93 

0-63 

45-28 

0-35 

966 

VEb 

54-62 

0-76 

44-35 

0*40 

10-28 

VEc 

54*88 

0-77 

4518 

0-43 

9-71 

VEd 

61-62 

1-34 

44-49 

0*79 

17-13 

VEe 

57-63 

118 

47-07 

0-98 

10-56 

8. Water ( WA) 

55-21 

0-44 

52*33 

0-40 

2*89 

WAa 

54-41 

0*27 

53-33 

0-28 

1 09 

WAb 

55*68 

0-42 

51 62 

0-35 

4-06 

WAc 

58*68 

1-77 

48-50 

1 54 

10-19 

9. Recreational (RC) 

63-03 

3-78 

45-61 

1-47 

17*42 

RCa 

61-40 

2 51 

44-66 

0-96 

16-74 

RCb 

63*16 

3*90 

44-95 

1-28 

18-21 

RCc 

71*79 

5*67 

4813 

2-31 

23-66 

RCd 

60-04 

3*43 

44-28 

1*08 

15-75 

10. Vacant (VA) 

60-88 

2-27 

44-35 

0*93 

16*53 


the only land cover type that shows the highest irradiance values at night 
(52-54 Wm 2 ), while agriculture exhibits the lowest at night. In other words, all the 
different categories of land cover are cooler at night than during the day. By virtue 
of their homogeneity, the land cover types of water, vegetation, and agriculture 
exhibit small SDs in their irradiance values at night as compared with commercial, 
services, and industrial land cover types. Residential uses, which are again inter- 
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mediate in position in irradiance at night (45-46 Win 2 ), show much less variation 
than that during the day. 

Because each broad class of land cover comprises sub-classes of distinctive 
characteristics, their irradiance values vary. In residential uses, multi-unit dwellings 
(such as apartments and condominiums, coded as REd in tables 2 and 3) show 
higher daytime irradiance value (69 Wm" 2 ) than single-family houses with tree covers 
(such as, REa with >50 per cent of tree canopy and REb with 50-10 per cent tree 
canopy, which exhibit 62 and 65 Wm -2 respectively). The higher percentage of tree 
cover results in lower irradiance during the daytime. At night, the difference in 
irradiance among the four subclasses of residential uses is much smaller than that 
in the daytime (varying from 45 to 46 Wm 2 only). In agricultural uses, harvested 
and fallowed land (AGb) shows much higher irradiance in the daytime than row 
crop land (AGa) and pasture land (AGc) (table 3). At night, while all show lower 
irradiance, the contrast between harvested fields and cropped areas is very much 
reduced ( 42-44 Wm 2 ). As for the vegetation cover, there is a predominance of 
mixed coniferous and deciduous trees (VEc) in the city, mostly on the mountains 
and along the rivers. They show higher irradiance (55Wm“ 2 ) in the daytime and 
much lower irradiance (45 Wm -2 ) at night. Grass cover (VEd) shows the highest 
irradiance (62 Wm 2 ) during the day but much lower irradiance (45 Wm -2 ) at night. 
As a result, grass cover registers the highest day-night irradiance differential 
(17Wm 2 ) in the vegetation land cover category. The most noteworthy is the 
wetland vegetation (VEe), which is intermediate in day-time irradiance between 
forests and grass, but it shows much higher irradiance than other types of vegetation 
at night probably because of its association with water. 

For services land cover types, schools and universities (SEc) show higher irradi- 
ance (70 Wm -2 ) during the day than other services classes (table 3). The government 
uses (SEd, mostly located in Marshall Space Flight Center and Redstone Arsenal) 
also exhibit high daytime irradiance (67 Wm -2 ). At night, the irradiance values of 
these uses are much more similar (46 Wm -2 ). It should be noted that offices housed 
in high-rise buildings (SEa) register lower irradiance values (65 Wm" 2 ) than offices 
housed in low buildings (SEb) (69 Wm 2 ), suggesting that landscaping may have 
helped to lower the irradiance. For commercial uses, large shopping malls (COa) 
show slightly higher irradiance (72 Wm -2 ) than small shopping centres (COb) 
(70 Wm' 2 ) in the daytime, but both types of commercial uses show much lower but 
very similar irradiance at night (46-47 Wm -2 ; table 3). Industrial uses exhibit higher 
irradiance for light industrial uses (MNa with 69 Wm -2 ) but lower irradiance for 
heavy industrial uses (MNb) and warehouse uses (MNc) (63 and 64Wm -2 respect- 
ively). At night, heavy industrial uses (MNa) show the highest irradiance of the three 
with 48WtrT 2 as compared with 43Wm -2 for warehouse uses and 45Wm -2 for 
light industrial uses, so that there are some significant night-time differences among 
these three classes of industrial use (table 3). 

Many of these commercial and industrial areas have buildings with aluminum 
roofs which possess very low emissivities (i.e., very high reflectance). Hence, they 
cool off very rapidly at night. In contrast with the services, commercial, and industrial 
land cover types, water (with emissivity close to unity) shows little difference in 
irradiance during the day and at night. However, wetland ( WAc) exhibits the highest 
irradiance (59 Wm -2 ) in the daytime and the lowest irradiance (49 Wm -2 ) at night. 
Lakes (WAb) show slightly higher irradiance (56 Wm" 2 ) than rivers and streams 
(WAa) (55 Wm -2 ) during the day but lower irradiance (52 Wm -2 ) than rivers and 
streams (53 Wm 2 ) at night. 
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5.2. Relations between NDVI and irradiance of land cover types 

The value of NDVI varies from — 1 to -hi indicating the amount of vegetation 
(biomass) that is found in the land cover polygon. In this research, the relation 
between the irradiance of the land cover type and the amount of vegetation (as 
represented by NDVI) is investigated. The Pearson’s correlation coefficients were 
computed between the NDVI of each broad land cover type and its irradiance. The 
significance of each correlation coefficient was determined using a one-tail Student’s 
f-test. Table 4 indicates that NDVI values tend to be correlated negatively with the 
irradiance of all broad land use types for both day and night with the exception of 
water, which is correlated positively with the irradiance during the day, but negatively 
at night. The highest negative correlation during the day is with the irradiance of 
transitional or vacant ( — 094), residential ( — 0 86), agricultural ( — 0-84), and vegeta- 
tion (—0*74) land, all at 99-9 per cent or higher level of confidence. Services uses 
also exhibit a lower but equally significant correlation (—0-75). At night, NDVI 
shows the strongest negative correlation with the irradiance of recreation land 
( — 0-91) at 99*9 per cent level of confidence, although significant (at 99 9 per cent 
level ) negative correlations are also observed for water (—0*69), agricultural (—0*60), 
vegetation ( — 0*58), and services (—0*51 ) uses. The fact that the irradiance values of 
services uses exhibit significant correlation with NDVI values for both day and night 
indicates the occurrence of a green environment for the location of offices in 
Huntsville. 

The NDVI’s very strong (but negatively correlated) relationship with the irradi- 
ance of residential, agricultural, vegetation, and vacant/transitional uses implies that 
the higher vegetation amount which characterize these uses can bring down the 
surface temperatures in that land cover type. In other words, residential, agricultural, 
vegetation, and some vacant/transitional land uses in the city of Huntsville are 
conducive to the lowering of daytime surface temperatures. At night, residential land 
is less effective than agricultural and transitional land in lowering the surface temper- 
ature as indicated by the insignificant negative correlation between residential use 
with NDVI at night (table 4). 

Plots of commercial, industrial, and transportation uses against NDVI (figure 3 ) 
are distinctly different from those for residential, agricultural, and vacant/transitional 

Table 4. Pearson's correlation coefficients between NDVI and day- and night-time irradiance 
of the broad classes of land cover in the city of Huntsville, Alabama, U.S.A. 


Land cover class 

Number of 
polygons 

R 

(daytime) 

R 

(night-time) 

1. Agricultural 

64 

— 0*8367* 

-0-6022* 

2. Commercial 

31 

-0-3556 

-01 460 

3. Industrial 

21 

-0-4674 

-0-1946 

4. Recreational 

11 

-0*6799 

-0-9073 

5. Residential 

68 

-0-8590* 

-0-2416 

6. Services 

41 

-0*7470* 

-0-5047* 

7. Transportation 

12 

-01 069 

-0-4850 

8. Vacant/transitional 

20 

-0*9435* 

-0-2031 

9. Vegetation 

57 

-0*7428* 

-0-5759* 

10. Water 

24 

+ 0*6756* 

-0-6943* 


* indicates significance level at 99-9 per cent determined by Student’s 
f-test. 
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Figure 3. Scatter plots of NDVl versus irradiance from (a) commercial; (b) industrial; and 

(c) transportation land cover types. 


uses (figure 4) in that they occupy the quadrant with negative NDVIs, an indication 
that they do not have a lot of vegetation amount to help lower their surface 
temperatures. This is particularly worse for the commercial land use. 

Services that show a moderately strong negative correlation with NDVI for both 
day and night are in a mixed group. In some types of services uses, such as offices 
in low-density areas (the research park) as well as schools and universities, some 
vegetation is associated with them, which gives rise to higher NDVI. The military 
and government uses included in this broad category of services are more mixed in 
character, which range from office buildings in Marshall Space Flight Center to 
missile testing ranges in Redstone Arsenal. But as a group, services uses are associated 
with much more vegetation than commercial and industrial uses. 

The very strong relationship between NDVI and residential and agricultural uses 
prompts the authors to perform a linear regression analysis. For the relationship 
between NDVI and residential uses, the adjusted R 2 value is 0-73 based on 68 
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Figure 4. Scatter plots of NDVI versus irradiance from («) agricultural; (/)) residential; and 

(c) vegetation land cover types. 


polygons. The linear model developed takes the following form: 


or conversely: 


1^ = 66-87 -22-30 (X ndvi ) 


( 2 ) 


y ndvi = 2-23-0-03(X lre ) (3) 

where 7 ire or X in and X ndvi or y ndvj are irradiance and NDVI for residential land 
cover respectively. 

Similarly, for agricultural uses, the adjusted R 2 value is 0-70 based on 64 polygons. 
The linear regression equation is: 

y; aq = 62-87 -13-91 (X ndvi ) 


(4) 
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or conversely, 

y ndvl = 3-26-005 (X iaq ) (5) 

where Y iag or X iag , and or T ndvi are irradiance and NDVI for agricultural land 
cover type respectively. 

By using these equations, the irradiance of residential and agricultural uses can 
be predicted based on NDVI (or vice versa). Regression equations such as these are 
useful because thermal infrared data are not always available, while NDVIs which 
are derived from image data in the visible and reflected infrared bands are more 
easy to obtain. Thus, NDVI is a surrogate of the irradiance of a land cover type, 
from which some ideas of the surface temperatures of the land cover type can be 
obtained. 

6. Development of a spatial model of an urban heat island for Huntsville using GIS 

A land use/cover map of the city of Huntsville can be obtained from a visual 
interpretation of the C1R aerial photographs (figure 5). This map was digitized and 
then converted into the raster format and georeferenced. Because each land cover 
type has a distinctive thermal signature during the day and at night, as revealed by 
the findings in this research, the digital land use/cover map can be converted into a 
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Figure 5. Land use/cover map of Huntsville. The UTM co-ordinates of the map extent are: 
(I) upper corner: 525300mE and 3856700mN; (2) upper right corner: 550800mE 
and 3856700m N; (3) lower left corner: 525300m E and 3822800m N; and (4) lower 
right corner: 550800m E and 3822800m N. The live black pixels are control points 
used to register this map with other layers of data. 
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map of irradiance, one for the day (figure 6) and one for the night (figure 7), using 
the ‘assign attribute 1 function of the GIS. The degree of detail mapped depends on 
the level of classification scheme used for land use/cover interpretation. By using 
high-resolution aerial photography at the scale of 1:15000, very detailed land 
use/cover map (Level III/IV) can be obtained by manual interpretation. For simpli- 
city, this research just makes use of the average thermal signatures of the 10 broad 
classes of land use/cover of Huntsville (table 3). By drawing profiles (such as the 
N-S and E-W ones), with the aid of the GIS software, one sees the occurrences of 
numerous ‘peaks’ from the centre of the city outward during the daytime and at 
night (figures 8 and 9). However, during the night time, the ‘peaks’ in the day have 
become ‘troughs’ (figures 8 and 9). The night and day irradiance difference map 
reveals that these troughs are in the centre of the city where most commercial and 
service uses occur (figure 10). Huntsville is characterized by the occurrence of agricul- 
tural land cover (30-9 per cent) and forests (35-9 per cent) in the periphery as well 
as inside the city. Residential uses (16*8 per cent) are also associated with tree cover. 
Commercial (2*0 per cent) and services (7*6 per cent) uses are small compared to the 
other uses. The spatial pattern of these different land use/cover classes can affect the 
development of urban heat islands. It seems that the layout of land use/cover for 
the city of Huntsville favours the development of urban heat island not only as a 
difference in temperature between the city centre and its periphery, but also among 
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Figure 6. Daytime irradiance map (in Wm' 2 ) of Huntsville. The E W and N-S cross section 
lines for figures 8 and 9 are indicated. The five black pixels are control points used to 
register this map with other layers of data. 
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Figure 7. Night-time irradiance map (in Wm" 2 ) of Huntsville. The E-W and N-S cross 
section lines for figures 8 and 9 are indicated. The five black pixels are control points 
used to register this map with other layers of data. 


agricultural, vegetation, and commercial uses in the whole city, thus giving rise to 
the emergence of many small urban heat islands throughout the city. Such a develop- 
ment is enhanced by the topographic characteristics of Huntsville, which may have 
also induced radiation fog to the warmer residential areas at the lower level at night 
as revealed by the flight line of channel 13 image along the eastern part of the city 
where mountains predominate. 

7. Conclusions 

There are some important findings from this research. High-resolution thermal 
infrared image data acquired from the ATLAS sensor system after correction for 
atmospheric transmittance and path radiance followed by blackbodies temperature 
calibration produce radiance values in 32-bit accuracy for the assessment of the 
urban heat island effect. Thermal signatures of different land cover types in the city 
for day and night help to shed light on their roles in contributing to the urban heat 
island phenomenon. It was found that commercial, industrial, and services uses 
exhibit the highest daytime irradiance while water, agriculture, and vegetation the 
lowest. Water, which shows very small day-night difference in irradiance, stands out 
unique as being the coolest in the daytime and the warmest at night because it 
possesses very high emissivity. Residential uses occupy an intermediate position 
because of their varying degrees of association with tree cover. A study of the 
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Figure 8. E-W profiles across — daytime and --- night-time irradiance maps in figures 6 

and 7. 



N-S Distance steps (1 step = 571.5 m) 

Figure 9. N-S profiles across — daytime and --- night-time irradiance maps in figures 6 

and 7. 
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Figure 10. Night-day irradiance difference map (in Wm 2 ) of Huntsville. 

relationship between NDVI (as a surrogate of vegetation amount ) and the irradiance 
of each category of land cover (as a surrogate of surface temperature) reveals the 
importance of vacant/transitional, residential, agricultural, and vegetation land cover 
types in contributing towards lowering their surface temperatures by virtue of their 
association with biomass. Based on their strong negative correlations, regression 
equations can be developed to predict the irradiance of the residential and agricultural 
land cover types from NDVIs. 

Based on a broad land use/cover map interpreted from aerial photographs, a 
spatial model of urban heat island for Huntsville has been developed by substituting 
in the thermal signatures associated with the land cover types. The model suggests 
that favourable conditions exist in Huntsville for the development of pockets of day- 
night temperature differences either between the city centre and the periphery, or 
among contrasting land cover types inside the city. The spatial layout of the land 
use/cover in a city has a great impact on the development of heat islands. Therefore, 
urban planning can be applied to avert or alleviate the effect of urban heat islands 
(McPherson 1994). 

Finally, the high-resolution thermal infrared image data permit detailed character- 
ization of the highly complex urban land cover types and hence a more realistic 
assessment of the urban heat island effect. 
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